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ABSTRACT 

Context. CAL 83 is a prototype of the class of Super Soft X-ray Sources (SXS). It is a binary consisting of a low mass secondary 
that is transferring mass onto a white dwarf primary and is the only known SXS surrounded by an ionisation nebula, made up of the 
interstellar medium (ISM) ionised by the source itself. We study this nebula using integral field spectroscopy. 

Aims. The study of ionised material can inform us about the source that is responsible for the ionisation, in a way that is complemen- 
tary to studying the source directly. Since CAL 83 is the only SXS known with an ionisation nebula, we have an opportunity to see if 
such studies are as useful for SXSs as they have been for other X-ray ionised nebulae. We can use these data to compare to models of 
how CAL 83 should ionise its surroundings, based on what we know about the source emission spectrum and the physical conditions 
of the surrounding ISM. 

Methods. With the VIMOS integral field spectrograph we obtained spectra over a 25 x 25 "field of view, encompassing one quarter of 
the nebula. Emission line maps - H i, He n ,14686, [O ill] AA4959, 5007, [N n] ^6548, 5683, and [S n] AA6716, 673 1 - are produced 
in order to study the morphology of the ionised gas. We include CAL 83 on diagrams of various diagnostic ion ratios to compare it to 
other X-ray ionised sources. Finally we computed some simple models of the ionised gas around CAL 83 and compare the predicted 
to the observed spectra. 

Results. CAL 83 appears to have a fairly standard ionisation nebula as far as the morphology goes: the edges where H is recombining 
are strong in the low stage ionisation lines and the central, clumpy regions are stronger in the higher stage ionisation lines. But the 
He ii emission is unusual in being confined to one side of CAL 83 rather than being homogeneously distributed as with the other ions. 
We model the CAL 83 nebula with cloudy using model parameters for SXSs found in the literature. The Hen emission does not 
fit in with model predictions; in fact none of the models is able to fit the observed spectrum very well. 

Conclusions. The spectral line images of the region surrounding CAL 83 are revealing and instructive. However, more work on mod- 
elling the spectrum of the ionised gas is necessary, and especially for the high-ionisation level emission from CAL 83. In particular, 
we wish to know if the He n emission and the other nebular lines are powered by the same ionising source. 

Key words. Galaxies: Magellanic Clouds - ISM: Individual objects (CAL 83) - X-rays: ISM 



peratures in the range of ~2 - 6 x 10 5 K (kT * 17-50 eV) 
( [Rappaport, Pi Stefano, & Smith|1994b ). 

Van den Heuvel et al. (1992) were the first to propose a 
model for SXSs: an accreting white dwarf of mass 0.7-1. 2M 
accompanied by a normal star of mass 1.5-2M©. According to 
their model the supersoft X-ray emission and high luminosities 
are a result of steady nuclear burning of hydrogen accreted onto 
the white dwarf. The mass transfer from the main-sequence star 
to the white dwarf occurs on a thermal timescale via Roche-lobe 
overflow, with rates of ~ 1 -4 x lO" 7 M yr _1 . The dominant opti- 
cal light source is the white dwarf, which, with its accretion disc, 
completely outshines the donor star. The nuclear burning on the 
surface of the white dwarf allows it to retain the accreted mass, 
making it a possible proge nitor of accretion-induc ed collapse 
(Ivan den Heuvel et al. 11992b or Type la supernovae dRappaportl 



1. Introduction 

Between 1979 and 1981, the Columbia Astrophysics Laboratory 
(CAL) carried out a systematic soft X-ray survey of the large 
Magellanic cloud (LMC) using the Einstein Observatory. This 
resulted in the d etection of 97 X-ray sources ( [Lon g, Helfan d, & 
|Grabelsky|1981| ). Some of the detected sources are characterised 
by an unusually soft spectrum in which little or no radiation at 
energies above -0.5 keV is detected. Source no. 83, being one 
of a small number of such sources, became known as CAL 83 
(also known as LHG83, RXJ0543. 7-6822, and 1^0543.8-6823) 
and it is now regarded as the prototype of the class of super- 
soft X-ray sources (SXSs). In 1990, the first X-ray all-sky sur- 
vey was performed with ROSAT and many more similar sources 
were found in the Galaxy, LMC, an d M3 1 ([Kahabka & van den 
Heuvel 1997| |Parmar et al.||1998] ). These sources are charac- 
terised by a luminosity of ~ 10 J/ " 38 ergs -1 and effective tem- 



Send offprint requests to: pieter.gruyters@gmail.com 



|et al.|1994b[|Di Stefano & Nelson| 19961 |Hachisu et al.|1996| ). 

CAL 83 has a time-averaged luminosity L x (0.15-4.5 keV) 
of 3.2 x 10 37 ergs -1 and is identified with a variable B~16.8 
star with a blue continuum and strong narrow (-305 kms -1 ), 
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variable Heii/14686 emission (Co wley et al.| [T998). It is, in 
fact, a mass-transferring binary, with a secondary mass of ~ 
0.5M o ( [Cowley et al.|1998|) and a white dwarf prim ary mass of 
1.3 ± 0.3M G jAlcock et al.|l997l|Lanz et al.|2005). The orbital 
period of the system is 1.04 days ( Smale et al.||1988| . While 
CAL 83 is one of a dozen luminous, extremely soft SXSs now 
known in our Galaxy and the Magellanic Clouds, so far it is 
the only one known to have an ionisation nebula. Its inhomoge- 
neous nebula has bright emission extended over 37", with fainter 
emission extending out beyond twice that, and a total mass of 
~ 150M o ( [Remillard, Rappaport, & Macri|1995| = RRM95); ob- 
viously the nebula is ionised interstellar medium (ISM), rather 
than being intimately related to CAL 83 itself. 

SXSs emit copious quantities of photons in the range 20- 
200 eV, hence the radiation can ionise any gas surrounding them, 
creating an ionisation nebula. We expect these ionisation nebu- 
lae to be distinct from classic H n regions where the ionisation 
is the result of photo-ionisation by massive O-type stars by the 
absorption of higher energy photons while the lower energy pho- 
tons escape. In comparison, in the case of the SXSs it is the other 
way around: the lower energy photons do the ionising while the 
higher energy photons escape (RRM95). Models for SXS nebu- 
lae by Rappapo rt et al.| ( |1994a| ) predict that these should be dis- 
tinct from other astrophysical nebulae, in particular, [O in] and 
He ii /14686 should be far brighter than in classic H n regions, and 
the radial gradients of these and other lines much more gradual. 
These models have not yet been tested on real data of SXSs, and 
CAL 83 presents a good opportunity to do so. 



Using the VLT VIMOS ( |Le Fevre et al.|2003) integral field 
unit (IFU) we have observed one field around CAL 83, obtain- 
ing spatial and spectral information about its surrounding ISM. 
In this paper we report on our findings, presenting emission line 
flux maps, spectra, line ratio maps, and a comparison of our re- 
sults to models of i onised ISM created with the cloudy code 
dFerland et al.|1998| ). 



2. Observations and data reduction 

2.1. Observations 

CAL 83 was observed at the VLT of ESO using the VIMOS in- 
strument in its IFU mode. We used the small field-of-view (FoV) 
mode, resulting in a coverage of 27" x 27", covered by 1600 
spatial pixels (spaxels) at a spatial sampling of 0.67"per spaxel. 
The data were taken on 2005 Oct. 28 and Dec. 1. The spatial 
sampling is contiguous with a dead space between spaxels of 



less th an 10 per cent of the spaxel separation ([Zanichelli et al. 



2005J. We used the high resolution blue grism (HRblue, -0.54A 
pixel -1 ), covering a spectral range of 4 150-6200 A and the high 
resolution orange grism (HRorange, ~0.62A pixel -1 ) covering 
5250-7400A. The observing log can be found in Table [I] We 
have 2 exposures per grism, and there is a slight offset in the 
pointing for some of the observations. For the flux calibration 
we made use of two standard star observations, one for each res- 
olution. As is standard ESO practice, these are not (necessarily) 
observed on the same night as the astronomical target. Table [T] 
also includes the details of the data from which we derived our 
sky spectrum (Sec. 2.2.2). 



2.2. Data reduction 

The data were reduced and flux calibrated using the ESO 
pipeline GasgancQ which allows the user to organise calibration 
files and run pipeline tasks. The reduction involves three main 
tasks: (1) vmbias which creates the master bias frame; (2) vmifu- 
calib which determines the spectral extraction mask, wavelength 
calibration and the relative fibre transmission corrections; and 
(3) vmifustandard for producing the flux response curves from a 
summed spectrophotometric standard star spectrum. Finally, the 
output of these tasks are fed into vmifuscience which extracts 
the bias- subtracted, wavelength- and flux-calibrated, and rela- 
tive fibre transmission-corrected science spectra. The FoV of the 
VIMOS IFU is split into four quadrants, and the light from each 
quadrant goes to a separate CCD, therefore the data processing 
is performed separately on each quadrant, creating four fully cal- 
ibrated frames per science exposure. 

Before doing the flux calibration it was necessary to step 
out of the pipeline to improve on the fibre transmission correc- 
tion. Looking at the (Gaussian-profile fit) flux of the sky line at 
5577 A in each spaxel for each quadrant separately, clear sys- 
tematic intensity variations of the order 20-30% were found: the 
"spectrum" of flux versus spaxel number showed the same pat- 
tern for all the frames of the same grism. For each frame and 
each quadrant of each grism, we measured the deviation of each 
spaxel with respect to the overall median of the individual quad- 
rants. This value was then divided out of the science frames, 
reducing the variations to a random scatter of 5-10%. 

Finally, it was necessary to apply multiplicative corrections 
to the spectra of the four quadrants so that the quadrant-mean 
5577 A flux was the same for all quadrants. These rescaling fac- 
tors are 1.3, 0.9, 0.8 and 0.9 respectively for quadrants 1 to 4, 
with scatter between the different observations. 



2.2.1. Combining data 

The final data cube for each grism was produced by averag- 
ing the two exposures of each. The slight pointing offsets be- 
tween the frames had to be corrected for before the frames 
could be combined. We used the IRAF routine drizzle ( Mut chler] 
& Fruchter||1997| , which was developed for Hubble Space 
Telescope (HST) images, to do this. For each wavelength of our 
spectral range, we extracted the 2D (flux) image, and spatially 
resampled the one frame to match the world coordinate system 
(WCS) of the other frame, for the blue and orange data sepa- 
rately. The offsets to apply were determined from the three stars 
that are in our FoV, by measuring their position to a precision 
of 1/2 spaxel in either direction using continuum images con- 
structed from our data cubes and comparing between the images. 
The offsets are the same for all wavelengths (there is no evidence 
for differential atmospheric refraction in our data). Offsets are of 
the order 1.5-2.5 spaxels. As a result of these offsets, the final 
maps are slightly smaller than the IFU FoV, being 25.5" x 25.5". 
For the subsequent analyses we placed both cubes on the same 
artificial WCS grid (but with the correct spaxel relative positions 
and sizes) - the slight offset in pointing between the blue and 
orange cubes is hence noticeable in the figures in this paper. 



2.2.2. Sky subtraction 

Sky subtraction can be a problem with the VIMOS IFU, as there 
are no sky-dedicated fibers. In our FoV there are no spaxels with 
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Table 1. The VIMOS IFU observing log. 



Name 


Observation ID 


uaie 


Grism 


Exp. time (s) 


KA 


uec 


Bluel 
Orange 1 
Orange2 
Blue2 


206671 
206675 
206676 
206677 


28/10/2005 
01/12/2005 
01/12/2005 
01/12/2005 


HR blue 
HR orange 
HR orange 

HR blue 


1 x 1320 
1 x 1320 
1 x 1320 
1 x 1320 


05 45 33.34 
05 45 33.35 
05 45 33.11 
05 45 33.11 


-68 22 19.2 
-68 22 19.4 
-68 22 18.0 
-68 22 18.0 


Standard 1 
Standard 2 
Sky 


200152644 
200152643 
200152515 


28/10/2005 
28/10/2005 
12/10/2005 


HR blue 
HR orange 
HR orange 


1 x 123 
1 x 154 
1 x 123 


00 23 48.97 
00 23 48.97 
23 20 42.56 


-27 51 43.1 
-27 51 43.1 
+05 09 47.9 



only telluric emission. Therefore we looked in the ESO archive 
for long observations of standard stars taken with our IFU set- 
up and close in time to our observations. These frames are listed 
in Table 1 . The FoV of these frames contain only one star, cen- 
tred on each quadrant. We reduce the data through the pipeline, 
including flux calibration. The spaxels clear of stellar emission 
were then extracted. Ideally we would have one sky spectrum 
per quadrant, or even per spaxel, as there are slight variations 
in the spectral PSF with position on the CCD. However, lack of 
signal-to-noise ratio meant we had to sum all the sky data into 
a single template sky spectrum. To subtract this template from 
our data it was necessary to scale the template sky spectrum. We 
measured the line fluxes from the template and from each spaxel 
in our CAL 83 frames. There is no variation, outside of noise, in 
the sky fluxes from CAL 83 with spaxel (just as one would ex- 
pect) and hence we scaled the template sky to the median of the 
whole CAL 83 field. The sky lines chosen for scaling were not 
the brightest - AA5511, 5890, 5896 (and secondarily 6300, which 
can also have some Galactic contribution) - as this resulted in an 
under- subtraction of all other sky line fluxes. This may be due to 
variation in the telluric spectrum with time. Instead we scaled to 
a dozen or so mid-bright lines. The results of our sky subtrac- 
tion are shown in Fig.[T] In the plot showing the reddest spectral 
range, S -shaped residuals can be seen: this is a consequence of 
variations to the spectral PSF with position on the CCD which 
we could not correct for. 

We carried out the sky subtraction only for the orange data. 
For the blue data the telluric spectrum is less of a problem as 
there are no sky lines that can interfere with our CAL 83 spec- 
tral lines; and we are not interested in the continuum level, only 
the emission line fluxes. The only potential contamination can 
be at He 1 5876 A and around 5200 A (where lines of [Ni] may 
be found). Careful inspection of our data (in cube form, in row- 
stacked spectrum form, and comparison to a blue sky spectrum) 
shows that there is no Hei (redshifted: 5881 A) and no [N i] 
(-5205 A) emission above telluric. 



2.3. Emission line fitting and map making 

Our final frames contain 1600 spectra each. To fit spectra of 
this quantity we adopted an IDL-base d semi- autom ated fitting 
procedure called PAN (Peak ANalysis; [Dimeo 2005 ), modified 
by |Westmoquette et al.| ( |2007] ) to accept row- stacked- spectrum- 
format FITS files. The user chooses a spectrum in the stack, a 
spectral region to fit, and defines the models - usually a poly- 
nomial to the continuum and one to a few Gaussian profiles to 
emission lines - and the software proceeds to fit all spectra in 
the input frame. The Levenberg-Marquardt technique is imple- 
mented to solve the least- squares minimisation problem during 
the fitting. 





Sky 
— CAL83 







Residual 



"5% 



5900 
Wavelength (A) 




6800 
Wavelength (A) 



6850 



6900 



Fig. 1. Sky subtraction. Shown is the spectrum of a single spaxel 
of CAL 83 with the scaled telluric spectrum overplotted. Below 
each plot are the residuals. On the top one notices the difference 
in the flux of the brightest sky lines with the rest (see text). On 
the bottom one can see S-shaped residuals which are a result of 
spectral PSF differences (see text). 



For PAN to produce its best result a realistic error spectrum 
must be attached to each source spectrum. These error spec- 
tra were created by us in the following way: the raw science 
frames (with signal values of ADU) were converted to (assumed 
Poissonian) percentage error frames, adjusted and combined in 
the same way that the science frames had been, and then multi- 
plied into the final CAL 83 frames to create error frames in phys- 
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ical units. These were attached to the CAL 83 FITS files and read 
into PAN. 

PAN produces residual spectra (calculated as r t = (yf lt - 
yf ata )/cr h where <x; are the data uncertainties), and the model 
results (for the Gaussians these are the central wavelength, inte- 
grated flux and FWHM, plus their errors), the spaxel numbers, 
and the x 2 values are output to a file. Inspection of the residual 
spectra, the errors and^ 2 values, allows one to discard doubtful 
results. 

These fits were later to be turned into emission line maps. To 
do this one requires only the table (provided by ESO) that relates 
the spaxel number to its relative sky coordinates. Constructing 
the maps, displaying as greyscale or contour, and plotting maps 
on top of each other (contour over greyscale and accounting 
for shifts in the WCS) were done from the PAN-output files 
with /ython scripts and using the tools provid ed by HI PE, the 
Herschel Interactive Processing Environment ( Ott|2010 ). 

Note that the data were not dereddened before being fit in 
PAN. 



2.4. Dereddening 

To correct for reddening we used the Galactic reddening law of 
Howarth (1983 ) with the c(Hj3) values derived from our observed 
Galactic Ha/Hj3 li ne ratios in conjunction with the theoretical 
Case B ratio from Osterbrock (1989 ), which is 2.85 for a tem- 
perature of 10, 000 K. We calculated a c(Hfi) value for all spaxels 
in the FoV using the Galactic Ha/Hj3 ratio. We then dereddened 
our observed LMC Ha and H/3 fluxes (the LMC redshift sep- 
arates them from the Galactic lines) and calculated the LMC 
c(H/3) values. This, however, resulted in a FoV with negative 
c(H(3) values. As the map of our Galactic c(H(3) values is flat, 
we opted to use one median Galactic c(H/3) value for the whole 
FoV to deredden the data (i.e. we assume all the reddening is 
Galactic). The median ratio for Ha/H/3 derived from all the good 
signal-to-noise spaxels in our FoV is 3. 8 + 0.5, correspondin g to 
c(H/3) = which by following |Kaler & Lutz|fl985| ) and 

using E(B - V) = c(H/3) * (0.61 + (0.024 * c(H/3))), is avalue of 
E(B-V) = 0.23. This value is however fairly high and cannot be 
reconciled with values found i n the literature, i.e., the E(B - V) 
value found by |Kovacs ( 2000| ) is 0.13,[Oestreicher & Schmidt- 
^a!er| ( [T996l ) deduced a range of E(B - V) of (0-0.15) for the 
Galactic foreground extinction toward the LMC, and RRM95 
adopted E(B-V) = 0. 1 in their paper about the nebula surround- 
ing CAL 83. However if we adopt the lower limit for the c(H/3) 
value we find a value of E(B - V) = 0.11, in perfect agreement 
with the literature. Hence we adopt c(H/3) = 0.18 to deredden 
our fluxes. We repeat, we chose this lower limit mainly to make 
comparison to previous CAL 83 data more straightforward, and 
because the lower limit is within the standard deviation (scatter) 
in our field of view. 

All values quoted henceforth (in tables, images and plots) are 
the dereddened values. The spectra presented in the paper are not 
corrected for extinction. 



3. Results 

With the aid of emission line maps we can discuss what the 
CAL 83 nebula looks like. We will compare these emission maps 
to previous observations of the nebula. We have also extracted 
the global spectrum of the nebula and of the CAL 83 stellar 
source, and will present the latter and briefly compare the stellar 
He ii line to what is already known of this feature. 



3.1. Global spectrum 

Fig. 2 shows our spectrum of the ionisation nebula surrounding 
CAL 83 taken from our 25.5" x 25.5" FoV, which corresponds to 
an actual size of 7.5 x 7.5 pc at the distance of the LMC (55kpc; 
|Smale et al.fT988] ). The stars (CAL 83 itself and two field stars) 
have been excluded from this spectrum. The emission lines in the 
spectrum are labelled. Of particular interest is the detection of 
He ii /14686 which has never been detected in the nebula before. 
In addition to He n /14686 we detected the usual forbidden and 
Balmer lines. We also detected [O i] /L16300, 6363 and Balmer 
emission from our own Galaxy. From this spectrum, we mea- 
sured the emission line fluxes using the ELF package in DIPSO 
( [Howarth et al.|[2004| ). The results of the Gaussian fits to the 
emission lines are given in Tables[2]and[3] There we quote the in- 
tegrated reddened and dereddened fluxes derived from all spax- 
els without stellar emission in them, i.e., only the ISM spectrum 
(summing 1213 and 1233 spaxels together respectively for the 
blue [4300-5100 A] and orange [6300-6800 A] emission lines). 
Along with the integrated fluxes we also present the average flux 
per 0.67 x 0.67 "spaxel. The quoted errors in the tables include 
the flux and fitting errors. Reddening errors are not included, as 
is standard practice. To give a feeling for the order of magnitude 
of the reddening errors, for H/3 and Ha, with reddening errors 
included the fluxes would be (387 ± 75) x 10" 16 erg cm -2 s" 1 and 
(968 ± 213) x 10" 16 erg cm" 2 s" 1 . 



3.2. Emission line morphology 

From the PAN fit results emission line maps were made. These 
are shown in Fig. [4] The two brightest lines are included - [O in] 
and Ha (the H/3 map looks much as the Ha) - and in Fig. [3] 
our maps are compared to the images of RRM95. We see a 
good correspondence between our and RRM95's observations; 
and see that we observed about 1/4 of the full nebula. The 
lower panels in Figj4] show the ions of lower ionisation stage 
- [O i] ,16300, and [S n] ,16716 ([N n] ,16583 shows similar be- 
haviour as [S ii] /16716, and hence is omitted). 

RRM95 discussed the morphology of the CAL 83 nebula as 
seen in their images. In both ions it is a bright, globally sym- 
metric but irregular nebula that extends outwards from CAL 83 
the stellar source. In [O m] the inner region is brightest in an arc 
which reaches a peak intensity near a radius of about 3.3 pc and 
extends to a radial distance of 7.5 pc from CAL 83. There is also 
fainter emission visible out to ~25 pc from CAL 83. The arc can 
also be seen in Ha, but it is less pronounced. The region inside 
the arc at < 1 pc from CAL 83 is a local minimum in [O m] 
and less so in Ha. Our FoV includes the inner and middle neb- 
ular regions, which are yellow and red on the RRM95 images. 
RRM95 find a ratio of 5:3 for the yellow:red for both mapped 
ions, which is close to what we find. On our maps CAL 83 is in 
the SE corner. 

The position of CAL 83, indicated by a white star on the 
maps, was derived by comparing a continuum map of our FoV 
with the Ha image of RRM95, who identified CAL 83 as the star 
in the centre of the nebula. Except for H i, none of the ions shown 
in Fig.|4] peaks at or around the position of CAL 83 the stellar 
source. This peak in H i at CAL 83 is due to emission from the 
star itself (as we show later). For the Ha map, uniquely among 
all ions, one can see that at the very SE corner of our FoV (left 
of the star) the emission is brightening (and we show later that it 
is here that the nebular He n emission can be found). It is unfor- 
tunate that our FoV is cut off here. Otherwise the morphology in 
Ha and [O m] is similar, as can be seen by comparing the colour 
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Fig. 2. Spectrum of the ionisation nebula of CAL 83 taken from our 25.5"x 25.5" (7.5 x 7.5 pc 2 ) FoV, excluding the contribution 
from the three stars (2 field stars and CAL 83). The spectra were extracted from 1213 and 1233 spaxels for the blue and red part of 
the spectrum, so the area from which the spectra are extracted differs between the blue and the orange observations. Note that the 
spectra are not dereddened. 



Table 2. Emission line measurements for the summed spectrum of the CAL 83 nebula. 



Line 




fW 


F(A) 


1(A) 


no. spaxels 


7 Wspaxel x 10 " 3 




(A) 








(#) 


Hy 


4340.5 


1.157 


81.8 + 9.8 


130.0+15.6 


1213 


107.1 + 12.9 


Hen 


4685.7 


1.043 


48.1 + 3.1 


73.8 + 4.8 


1213 


60.8 + 4.0 


H/3 


4861.3 


1.000 


257.0 + 3.7 


387.6 + 5.6 


1213 


319.5 + 4.7 


[O in] 


4958.9 


0.976 


224.5 + 6.8 


335.2 + 10.2 


1213 


276.3 + 8.4 


[O m] 


5006.8 


0.964 


652.6 + 8.0 


969.8 + 11.9 


1213 


799.5 + 9.8 


[Oi] 


6300.3 


0.718 


152.1 + 13.4 


204.3 + 18.1 


1233 


165.4 + 11.4 


[Nil] 


6548.1 


0.682 


60.9 + 10.6 


80.6 + 14.0 


1233 


65.4 + 11.4 


Ha 


6562.8 


0.680 


609.2 + 10.6 


968.2 + 14.0 


1233 


785.3 + 11.3 


[Nil] 


6583.5 


0.677 


111.4+10.6 


147.1 + 14.0 


1233 


119.3 + 11.3 


[S II] 


6716.4 


0.658 


287.2 + 9.9 


376.4 + 13.0 


1233 


305.3 + 10.5 


[S n] 


6732.8 


0.656 


194.4 + 9.8 


254.6 + 12.9 


1233 


206.5 + 10.4 



Notes. Absolute observed reddened (F(A)) and dereddened fluxes (1(A)) are given with their corresponding errors, in units of 10 16 ergcm 2 s 1 . 
The adopted reddening curve, f(A) and the average flux per spaxel ^/0 S p aX el are a ^ so gi yen (0.67 // ). The value of c(Hfi) is 0.18. 



image to the overplotted contours: a bright peak in the centre of 
our FoV and secondary peaks along the E and S edges. The peak 
in [O m] is located closer to CAL 83 - the ionising source - than 
that in Ha. 

The other ions shown - [O i] and [S n] - also peak in the 
same general area (the centre of our FoV) but shifted slightly 
further away from CAL 83 than Ha, and they are faint in the 



radial region around CAL 83 itself (i.e., in the SE quadrant). In 
[S n] the morphology almost appears to be that of (part of) a ring 
of emission with CAL 83 at the centre. 

The order of the flux peaks in our nebula, with respect to 
the distance from CAL 83 [the ionising source], is: Hen, [O m], 
Ha, and then [N n], [S n] (and [O i]) peaking farthest out. This 
is in the order of their ionisation potentials, being 54.42, 35.12, 
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Table 3. Possible emission lines measured from the summed spectrum of the CAL 83 nebula. 



&obs 




F(A) 


1(A) 


no. spaxels 


^spaxel x 10 


(A) 








(#) 




A A C\C\ A AC A* 

44UU— 44j4 


1 1 rv^ 


4j /.j 


/zU.o 


1213 


CO/l 1 

jb>4.1 


4659.4* 


1.050 


61.9 + 17.4 


95.3 + 26.8 


1213 


78.6 + 22.1 


4663.9* 


1.048 


14.8 + 15.2 


22.8 + 23.4 


1213 


18.8 + 19.3 


6462.6 


0.694 


70.6 + 10.6 


93.9 + 14.0 


1233 


76.2 + 11.4 


6495.3 


0.689 


55.5 + 10.6 


73.7 + 14.1 


1233 


59.7 + 11.4 


6532.7 


0.684 


38.5 + 10.6 


51.0+14.0 


1233 


41.4 + 11.3 


6634.7 


0.670 


61.7 + 10.6 


81.3 + 13.9 


1233 


65.9 + 11.3 


6683.8 


0.663 


29.1 + 8.2 


38.3 + 10.7 


1233 


31.0 + 8.7 



Notes. Absolute observed reddened (F(A)) and dereddened fluxes (1(A)) are given with their corresponding errors, in units of 10 16 ergcm 2 s 1 . 
The adopted redding curve, f(A) and the average flux per spaxel ^(^) S p aX el are a ^ so gi yen (0.67"). Lines marked with * are blends. The value of 
c(H0)isO.18. 




13.60 and 14.54/10.36, respectively. Looking at the overall mor- 
phology in our 1/4 field of view of the nebula, we see that the 
[N n], [O i] and [S n] zones are essentially coincident with each 
other (with [S n] being the brightest), and while they also overlap 
with the regions of the brightest Ha and [O in] emission, they are 
more-or-less absent in the immediate area around CAL 83. So far 
this nebula has the appearance of a more-or-less standard ionisa- 
tion nebula with a single ionising source of high temperature. 

Note: the slight shift in FoV for the ions from the blue and 
the orange maps is due to the slight offset in the pointing for 
these observations. Noisy edge spaxels have been blanked out. 

3.3. Hen nebular morphology 

The most interesting thing about our observations is the detec- 
tion of He H/14686 in the nebula of CAL 83, as it has never 
been observed before. The Heii/14686 emission line has been 
detected in the stellar spectrum (see Sec. |3.4| ) and it has been 
predicted that it should be found in the nebulae of SXSs (e.g., 
Rappapo rt et al.| 1994a ). The stellar line is discussed in the next 
section, here we discuss the He n morphology. 

The map of He n is presented in Fig. [5] We could only mea- 
sure this emission line in each spaxel for one quadrant, the one 
containing CAL 83 itself. On the map one can see a bright re- 
gion, which is where CAL 83 is located. To confirm this, and to 
check whether this bright emission arises from a point source 
(i.e., CAL 83 itself) we compared the spatial distribution of He n 



at the location of CAL 83 to the spatial distribution of CAL 83 
from a continuum spectral region (i.e., from the star itself). 
These continuum contours are shown also in Fig. [5] It is clear 
that the brightest emission in He n arises from a point source. 

There is faint, but obvious Hen emission to the immedi- 
ate south of CAL 83 (where Ha also shows a peak; see Fig. [4]). 
Although we have not, unfortunately, mapped the entire sur- 
roundings of CAL 83, it is clear that this Hen emission is not 
symmetrically distributed about CAL 83 and is fairly confined. 
With an instrumental broadening of ~ 100 kms -1 , the nebular 
line is spectrally unresolved and shows no radial velocity gradi- 
ent that we can measure. The flux of He n coming exclusively 
from the bright Hen nebula (see L-shape in Fig.[7]), excluding 
CAL 83 itself, is (29.98 ± 1 .4) x 10" 16 erg cm" 2 s" 1 . Even though 
this is the brightest zone in the He n nebula, this is not all the 
He H/14686 flux in the FoV-the total flux measured over our 
FoV is (73.8 ± 4.8) x 1 0" 16 erg cm" 2 s" 1 -but unfortunately we 
could not measure the line per spaxel outside of the area shown 
in Fig.[5j 

We note that we did not detect any He i in our FoV. 



3.4. Hen (and Ha) stellar line 

We have extracted the stellar spectrum of CAL 83 from the six 
spaxels that contain the flux from the central star, CAL 83. The 
ISM contribution has been removed from this full-wavelength 
range spectrum, shown on the left of Fig. [6] In the right of Fig.|6] 



6 



P. Gruyters: The ionised nebula surrounding CAL 83 




Fig. 4. Top: [Om] map (left) and Ha map with overplotted [Om] contours (right). The colour ranges are 0.17-0.98 for Ha and 
0.19-1.85 for [Om], in units of 10 _16 ergcnr 2 s _1 . Bottom: [Oi] ,16300 (left), and [Sn] ,16716 (right) maps, all with Ha contours 
overplotted. In the [S n] image we have indicated the transition region over which the [S n] over [N n] ratio increases from 2.3 to 



3.5 see section 



4.2 



The colour ranges are 0.015-0.288 for [Oi] and 0.106-0.454 for [S n], in units of 10 



" 16 erg cm 2 s 1 . Colour bar 



at the bottom indicates the normalised colour range. E is top and Nis right, and the FoV is 25.5" x 25.5" (7.5 x 7.5 pc 2 ). The spatial 
PSF is about 3 spaxels in height and 2 spaxels width. 



we show a zoom-in on Hen and Ha, which are the only two 
emission lines in our spectrum that are also bright peaks at 
the position of the star. For this zoom, the ISM continuum has 
not been removed: instead we show for comparison the area- 
normalised spectrum of the ISM in these two emission lines. The 
stellar line is much broader than the nebular line and has more 
emission in its blue wing. While we cannot unambiguously es- 
tablish whether there is any nebular emission at the location of 
the stellar source, it is clear that any such contamination is small 
for He ii but more important for Ha. Our measured fluxes are not 
corrected for any ISM contamination, as we cannot be sure that 
ISM emission is blended with stellar emission. 

The stellar spectrum i n the left panel o f Fig. [6] is si milar to 
the spectra published by Crampto n et al.| ([1987 ) and |Cowley| 
|et al.| ( [T998] ). On the blue side of He n we fi nd a blend of carbon 
lines (Cm, Civ), these being first noted by Pakul l et al. ([1985). 
This emission, together with Ovi at 5290 A, implies a very 
high degree of ionisation (IPqvi ~ eV). We also note the 
detection of the [Fe x] A6315 coronal line in emission as first 
observed by |Qliveira & Steiner| ([2006). 



Noting the variable position and strength of He n, Crampton 



et al. (1987) argued that Hen and the C and O emission fea- 



tures arise from the accretion disc of the CAL 83 system. This 
was later confirmed by Cowl eyet al.| ([T998), who finds evidence 
that the emission is linked to matter being ejected in bipolar 
flows, which may arise from a 69-day precessing accretion disk. 
For our observations we find FWHM ~305 km s" 1 agreeing well 
within the range of values found by |Smale et al.| ( [T9§8] ), i.e., 233 
to 376 km s" 1 over a time period of one year. The FWZM extends 
to a velocity of 1150kms _1 while Cowley et al.| ( |l998] ) state a 
velocity of ~ 2450 km s" 1 derived from a single observation. 



Smale etal. ( 1988| ) deduced a time-averaged luminosity from 
the Hen stellar line, L*(Hen) = 3.6 x 10 33 (d/55 kpc) 2 erg s" 1 . 
They observe a variation in the strength of He n of a factor of 4 
in equivalent width over a 1 year period. We measure a flux of 
(11.3 ± 0.2) x 10" 16 erg cm -2 s" 1 , and converting this to a lumi- 
nosity (L = 47id(cm) 2 F(erg cm -2 s" 1 ) and abopting a distance of 
55 kpc) we obtain L*(He n)= (4. 1 ± 0. 1) x 10 32 erg s" 1 . Although 
this is a (slight) lower limit to the actual stellar flux, our FoV 
just clips the edge of the star, it is about a factor of 10 less than 
the time averaged value. (Remember that the measured flux con- 
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4.2. Ion-Ion maps: morphology 



Fig. 5. He ii /14686 flux map, Gaussian-smoothed over a 3 x 3 
pixel box. Overlaid in black are the contours representing the 
FHWM and FWZM [deduced from the CAL 83 stellar contin- 
uum emission]. Note that we show here only one quadrant of our 
FoV, that where He n emission was detected. The colour ranges 
from 0.003 to 0.49 x 10 ' 
sents the entire normalised colour range. 



16 erg cm 2 s 1 . The colour bar repre- 



tains contributions from the star as well as from any ISM located 
between CAL 83 and us.) 

We also find an extended blue wing in Ha (FWHM- 
270 km s" 1 ). We measure a flux of (6.5+0.4) x 10" 16 erg cm" 2 s" 1 
from the Ha stellar line, corresponding to a luminosity of 
L*(Ha)= (2.4 ± 0.4) x 10 32 erg s" 1 . We could not find any pub- 
lished values to compare to. 



4. Nebular diagnostics 

Using certain emission line ratios we have placed limits on the 
electron temperature and density of the CAL 83 nebula. We then 
use various common emission line ratios to study the variation 
of ionisation across the nebula. 



4.1. T Q and n e 

For photoioinised nebulae one can use the ratio of the sulphur 
lines ([S n] A6116/[S n] A6131) to determine a density. The me- 
dian flux ratio over our FoV is 1.44 with an uncertainty of ±0.14 
and a standard deviation of 0.2: so at best any variation in the ra- 
tio over the Fo V is small, (an d we can also see that it is random). 
Following Osterbrock (1989), with an assumed electron temper- 
ature of 10,000 K, this corresponds to an electron density of less 
than 10 particles per cm 3 . The average ISM density in the inner 
nebula (i.e., characterised by a radius of ~7.5 pc) was estimated 
by RRM95 to be in range of 4-10 cm -3 , in agreement with our 
findings. We assumed an electron temperature as we do not have 
all the [O in] or [N n] lines which are necessary for a more self- 
consistent measurement, although the lack of [Om] /14363 (i.e., 
the 1-cr limit) leads to an upper limit T e - 14, 500 K. 



One can characterize the spectral energy distribution (SED) of 
a source that is ionising a nebula by comparing the emission 
line fluxes from a modelled nebula (with input parameters of 
the nebular temperature and density, as well as source SED) to 
what is observed. However, not only must the nebular condi- 
tions be known, but one needs to be sure that all ions stages 
have been measured: this is not always possible, especially when 
only a limited wavelength coverage is available. Alternatively, a 
qualitative measure of the hardness of the radiation that is ionis- 
ing the nebula can be obtained by comparing the emission from 
lines of ions the lower ionisation potentials (e.g. [O i], [N n] and 
[S n]), to the lines of ions with higher ionisation potential, such 
as He ii and [O in] . The diagnostic ratios we will use to do this are 
[O m] ^5007/HyS, [N ii] /16583/Ha and [S n] \ M6716, 673 1/Har. 
The first ratio is a good reddening-free indicator of the mean 
level of ionisation (radiation field strength), since [Om] arises 
from an ion of a relatively high ionisation potential (35.1eV for 
+ — > 2+ ). The other two ratios are related to gas ionised at a 
lower energy level: the [Sn] and [Nn] lines arise from ions of 
lower energies (S° -> S + =10.36eV, N° -> N + =14.54eV). 

Now we describe what we see in Figs|4]and|7j and relate this 
to the ionisation potentials of the ions. 

[S ii] will rise in flux when hv reaches 10.36 eV. Roughly 
speaking, one would expect [S n] to remain ionised throughout 
the transition region between H + and H° (IP 13.6ev), and this 
we do see: [S n] is found in the region where Ha lies but extends 
further out, where the Ha flux drops. 

Nitrogen has a first ionisation potential of 14.54 eV, slightly 
higher than that for Ha, and so, assuming that a transition region 
exists between singly ionised hydrogen and neutral hydrogen (as 
it does for CAL 83), the [S n] emission will become dominant 
over the [Nn] emission as one moves through this transition 
region going away from the ionising source ( [Saraph & Seaton 
|1970|) . We have some evidence of this: on the map of [S n] in 
Fig.|4J we have drawn an arc where, as one moves away from 
the ionising source, the mean [Sii]/[Nii] ratio increases from 
2.5 to 3.5. This border is coincident with the region where the 
Ha emission drops. 

To have [O i] emission one requires neutral oxygen, and photons 
of hv > 13.6eV (IP for ionising to the first level) will decrease 
its abundance and hence the [Oi] flux. Hence the [Oi] flux 
should be strongest farthest away from the ionising source: 
while [O i] is faint on our maps, it is found throughout the region 
where [S n] is, and it is likely that for [O i] (as well as [S n] 
and [Nil]) we have not observed the entirely of their emission 
regions, which will extend much further beyond the edge of our 
maps. 

The top left panel in Fig. [7] shows that the peak of the 
[O m]/Hj3 ratio lies near the central region that is bright in all 
ions, and that it is fainter around CAL 83 itself. The more highly 
ionised region, characterised by a high [O m]/Hj3 ratio, fills the 
"inner nebula" of RRM95 (i.e., the S-E diagonal half of our 
FoV). The regions of low ionisation states, characterised by rel- 
atively high [S n]/Ha and [N n]/Ha ratios (see Fig.[7]) are found 
over and beyond the edge heading towards the outer nebular re- 
gion of RRM95 (i.e. the N-W diagonal half of our FoV). The 
map of [O i] /16300/Ha, and also [S n][Ha, shows clearly that 
these ratios are greatest where the Ha flux is lowest: arising 
from the edges of the "nebula" where the ionisation of H is drop- 
ping. In objects which are photoionised by a spectrum contain- 
ing a large fraction of high-energy photons these zones in which 
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Wavelength (A) Velocity (km/s) 

Fig. 6. Left: CAL 83 stellar spectrum extracted from 6 spaxels that include the source. The ISM contribution has been removed 
and the detected emission lines are labelled. Right: a zoom on the stellar He n and Ha line. The spectrum is plotted in velocity 
space, normalised to the velocity of the LMC: 270 km/s. Shown for comparison are the ISM lines, normalised to the same area over 
which the stellar source extends, and taken from the area where He n emission can actually be found (see Fig. [5}. Note that the ISM 
emission has not been subtracted from the spectrum. The spectra are not dereddened. 



the ionisation of H is dropping are quite extended (Veille ux~&| 
|Osterbrock|1987] ). 



4.3. Ion-ion ratios: diagnostic diagrams 

We now take a more quantitative look at the line-ratio maps. 
Emission line-ratio values for interesting regions in the nebula 
can be found in Table [4j The regions are drawn on the up- 
per right plot in Fig. [7] ([S n]/Hor). The triangle- shaped region 
corresponds to the "faint region", the oval to the "central re- 
gion" and the L- shaped region to the "Hen region". We can 
compare our derived ratios with those of different astronomical 
sources by placing them in a diagnostic diagram. A diagnostic 
diagram helps to classify the dominant energy source in emis- 
sion line sources using optical spectral diagnostics in a qualita- 
tive way. The diagrams are based on empirical findings and can 
be compared quantitatively to theoretical expectations. Our di- 
agnostic diagrams combine a measure for the hardness of the 
ionising radiation ([O in] /15007/H/3) with measures for the low 
ionisation degree ([O i] /16300/Hor, [S n] M6716, 673 1/Hor and 
[N ii] ,16583/Ha) in the nebula. 

These diagnostic diagrams were designed to distinguish 
narrow-line active galactic nuclei (AGN) from Hn region-like 
objects but can be used for any gaseous emission line-objects. 
Rappapo rt et al.| ( [T994a| ) suggested another potentially useful di- 
agnostic diagram for SXSs, in particular to intercompare them 
among Hn region-like objects. This diagram is based on two 
important emission lines which should be found in the spectrum 
of an SXS: [O m] ,15007 and He n ,14686. 

We compare our results with the ratios found in ultralu- 
minous X-ray sources ( ULXs). The data for these ULX s were 
taken from the work of |Abolmasov et al.| ( [2007] [2008] ). ULXs 
can be regarded as the big brothers of SXSs, in the sense that 
they are harder X-ray sources, but radiate extreme luminosities 
of 10 39 " 41 erg s" 1 in the 0.5-10 keV band, with significant soft 



limit for a 10M© black hole, leading to the speculation that they 



may harbour intermediate-mass black ho les (1MB Hs ; [Colbert & 
|Mushotzky|1999[ [Miller & Colbert|2004[ ). However, subsequent 
analyses support super-Eddington emission scenarios f or stellar 
mass black holes in all but the brightest ULXs (e.g., IRoberts 



[2007] [Gladstone, Roberts, & Done1|2009] ). A number of ULXs 
in nearby galaxies have now been identified with large nebulae 
(e.g. [Pakull & Mirioni|[2002l jRoberts et aT1[2003l [Abolmasov 



comp onents in their X-ray spectra (e.g., Gladsto ne & Roberts 
2009 ). At these high luminosities they exceed the Eddington 



|et al. 12007) , some of which display characteristics of X-ray ion- 
isation in their central regions. They therefore present an inter- 
esting comparator for the nebula around CAL 83. 

We have plotted our results for different areas of the CAL 83 
nebula on the diagrams of Fig. [5] In each we have plotted values 
for the three regions marked in Fig. [7] ([S n]/Har) and given in 
Table [4| and one point reflecting the average value over our FoV 
(excluding the stars). The Hen region in the south-east corner 
of our FoV is interesting because that is where the only nebular 
He ii is found, as well as a peak in Ha; the faint region in the 
north-west corner corresponds to the edge of the inner nebula, 
where the material is dropping in ionisation level; and the cen- 
tral region of our FoV is where [O in] and Ha are bright. Note 
that CAL 83 is the first and only SXS around which an optical 
ionisation nebula has been detected, so this is the first time an 
SXS nebula has been placed on these diagnostic plots. 

The CAL 83 nebula shows characteristics of Hn region- 
like objects and AGNe. The [N n] 26583/Ha ratio charac- 
terises our nebula as an Hn region-like object, while the 
[S n] /U6716, 6731/Ha places it on the side of the AGN-like ob- 
jects suggesting it has a more extended zone of partially ionised 
hydrogen than Hn region-like objects since AGNe are charac- 
terised by a harder radiation field. The placement of the different 
regions of CAL 83 shows some spread, but the points always 
lie in the same general area of the plots as each other. There 
seems to be no great distinction between this SXS and the ULXs 
which is interesting since nebulae surrounding ULXs are created 
by hard X-ray radiation. In the He n A4686/HJ3 diagram in Fig.[8] 
we include theoretical values we calculate for model SXS nebu- 
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D 0.2 0.3 OA 0.5 0.B 0.7 O.H 0.B 1 



Fig. 7. Top left: Ratio map of [O in] A5007/HJ3 with the [O in] /15007 contours overplotted as reference and the position of CAL 83 
is noted with a white star. The colour range is 0.05-6.5. Top right: Ratio map of [S n] AA6116, 6731/Ha with the Ha contours 
overplotted. The colour range is 0.15-1.2. The white shapes show the regions from which we measured line ratios: the triangle 
corresponds to the faint nebula region, the oval to the central region and the L-shape to the He n region. Bottom left: Ratio map 
of [N n] /16583/Ha with the Ha contour overplotted. The colour range is 0.05-0.3. Bottom right: Ratio map of [O i] A6300/Ha 
overplotted by the Ha contours. The colour range is 0.03-0.47. The colour bar represents the normalised colour range. All maps are 
Gaussian-smoothed over a 3 x 3 pixel box. E is up and N is right and the FoV is 25.5" x 25.5" (7.5 x 7.5 pc 2 ). 



Table 4. Diagnostic emission line ratios 



Region (no. spaxels) 






i [S ii] 


i [N II] 


los^ 


He ii 

(xl(T 16 erg cm -2 s _1 ) 


w 

(xl0" 16 ergcnr 2 s" 1 ) 


Nebula (Full FoV: 1213/1233) 


0.39 


-0.68 


-0.19 


-0.82 


-0.72 


73.79 


387.6 


Hen region (SE: 79) 


0.49 


-0.82 


-0.31 


-0.85 


-0.23 


29.98 


55.78 


Central region (171) 


0.57 


-0.60 


-0.14 


-0.77 


-0.70 


11.66 


59.43 


Faint region (NW: 556) 


0.31 


-0.57 


-0.12 


-0.76 


-0.89 


18.43 


141.3 



Notes. The diagnostic ratios for different regions in the nebula along with the measured He n and HJ3 flux and the number of spaxels over which 
the fluxes were measured. The 'nebula' values correspond to the fluxes measured from our full FoV using 1213 and 1233 spaxels for respectively 
the blue and orange observations. 



lae (Sec. 5). The models are characterised by high values for the 
[O m] /15007/H/3 ratio and substantial He n /14686 emission. Our 
value from the He n region is consistent with the He n A4686/H/3 
range suggested by the models but we did not observe the pre- 
dicted high [O m] A5007/HJ3 flux ratios. 



5. Theoretical Models vs. Observations 



Rappaport et al. (1994a) created some simple models for SXS 



nebula, based on a spherical cloud ionised by a hot blackbody. 
The total abundances and abundance gradients were compared 
to the little data they had for CAL 83. However, their plots and 
tables were taken from models run with Galactic abundances. 
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Fig. 8. Diagnostic line-ratio diagrams for different regions in the CAL 83 nebula compared with ultra-luminous X-ray sources 
(ULX), which have much ha rder X-ray spectra, as well as star forming galaxies (SFG) and AGNe. ULX data were taken from 
Abolmasov et al.| ( [2007||2008] ). The name for each ULX is given next to the symbols. The positions of different parts of the CAL 83 
nebula are given by filled upright triangles. Solid curves represent the theoretical 'maximum starburst' line from Kew ley et al.| ( |200l| 
which distinguishes AGNe from Hn region-like objects (such as SFG). Left: [O ni] /15007/H/3 vs [S 11] AA6116, 6731/Ha. Right: 
[O m] A5007/HJ3 vs [N 11] ^6583/Ha. Bottom: [O in] /15OO7/H0 vs He 11 /14686/H£. The filled dots here correspond to our cloudy 
models for SXSs. The numbers near the filled circles indicate the model number from Table[5] For comparison we have included 
starburst galaxies and H 11 regions, the data for these objects being taken from the same paper as the SXS models. Arrows represent 
upper limits; vertical bars with horizontal arrows drawn adjacent to certain starburst galaxies and H 11 regions indicate common sets 
of upper limits. 



Since CAL 83 is in the LMC, we decided to run new models with 
LMC abundances. For this we use the latest version (c 10.00) of 
spectral synthesis code cloudy ( [Ferland et al.|19 98). In this sec- 
tion we will compare these models with our observations and 
discuss the results. Following Rappapo rt et al.| ( [l994a| ), we have 
produced 5 simple models for SXSs in which we differentiate be- 
tween luminosity, temperature and density. For comparison we 
have also included a model f or a classical H 11 region using the 
physical parameters given in Dopita & Evans (1986). The input 
parameters for the different models can be found in Table [5] 

The ionising source in all models is represented by a point 
source radiating as a simple blackbody with a given temperature. 
The point source is placed in the centre of a spherical symmet- 
ric and homogeneous gas cloud that has comparable elemental 
abundances relative to hydrogen of the LMC and were take from 



Table 5. Supersoft X-ray Source Model Parameters. 



Model 


r eff 


L 


n 


(#) 


(K) 


(ergs -1 ) 


(cm" 3 ) 


1 (standard) 


4x 10 5 


1 x 10 38 


10 


2 


2x 10 5 


1 x 10 38 


10 


3 


7x 10 5 


1 x 10 38 


10 


4 


4x 10 5 


1 x 10 37 


10 


5 


4x 10 5 


1 x 10 38 


1 


6 (Hn R) 


4.6 x 10 4 


5.1 x 10 38 


10 



|Rolleston et al.| ( |2002| ). The standard model (model 1) has an 
X-ray source of temperature 4 x 10 5 K (kT = 34 eV) and a lumi- 
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nosity of 1 x 10 38 erg s 1 , representative of cited values for SXSs 
in the literature ([Trumper et~aT 1991| Orio & Ogelman 1993 



Kahabka & van den I 



euve 



Hertz, Grindlay, & Bailyn||1993 ^ 
1997J). The hydrogen number density is taken to be 10 hydrogen 
atoms cm -3 . The other models are calculated by holding two 
of the parameters fixed with respect to the standard model and 
varying the third parameter. 



5.1. Emission line fluxes 

The output of cloudy includes intrinsic line intensities and lu- 
minosities of the cloud. The intrinsic emission includes all pro- 
cesses that affect the line formation and transfer. This includes 
collisional processes, fluorescence, line destruction by back- 
ground opacities such as dust or the Lyman continuum of hy- 
drogen, and recombination. The intrinsic intensities do not in- 
clude the effects of absorbers or scatters that do not lie within 
the line-formation region. 

Table [6] gives the values of the emission line fluxes for each 
of the six model nebulae and our measured values. All emission 
line intensities in the table are normalised to the corresponding 
HJ3 flux to correct for the fact that we only observed about l/4th 
of the full nebula. 

Our observed fluxes are those of Table 3, with l-cr upper 
limits for [Ni]/U5 199, 5202, [Oi] ,16363 and He i ,16678, as 
measured from the full FoV spectra using the ELF package in 
DIPSO dHowarth et aL|2004) . 



Unfortunately none of the models is a particularly good 
match to our observations. For most of the ions the models pre- 
dict larger line ratios (compared with Hfi) than we observe. The 
largest discrepancies are: (i) the extremely high [O iii]/H/3 flux 
ratio is not observed, where the measured values are typically ~ 
1/4 the SXS model predictions and about twice the Hn region 
model predictions; (ii) our observed Hen/H/3 ratio for the whole 
nebula is less than half the model predictions, note however that 
if we just take the Hen region (Hen/H/3 = 0.59), the ratio does 
match up; (iii) our observed [N n] ratios are 5 times to about an 
order of magnitude lower than the model predictions, and (iv) 
the [S n] ratios are a factor of ~3 lower than predicted by our 
standard model. 

We do not have a ready explanation for these systematic dis- 
crepancies. It is possible that the [O m]/Hf3 ratio in particular is 
affected by the fact that most of the [O in] (as demonstrated by 
RRM95) is outside our FoV. Indeed, the spatial inhomogeneity 
of the CAL 83 nebula, and the fact that our observations cover 
only a fraction of the structure, must have implications for our 
comparison of the simple model to the data that could only be 
resolved by a spatially complete analysis of the nebula. A more 
speculative explanation for the discrepancies could be that the 
CAL 83 nebula has unusual abundances for the LMC. This could 
work toward explaining the low nitrogen ratio we observe com- 
pared to the sulphur ratio. It could also be that we still do not 
have a full understanding of the astrophysics of these regions. 
Further observations of different parts of the nebula will provide 
us with better constraints on our models and may eventually aid 
us in better understanding the physics of these unusual nebulae. 



6. Conclusions 

In this paper, we have presented flux maps made from the fit- 
ting the emission lines of [O m] ,15007, [O i] ,16300, Ha and 



[S n] /16716 showing the morphology of these ions in our FoV. 
The morphology in [O iii] ,15007 and Ha match those found by 
RRM95 who imaged the full nebula in these ions. We find the 
edge of the RRM95 nebula is where low ionisation ions peak, 
while inside this the ionisation state is higher. 

We have also presented, for the first time, evidence of an 
He H/14686 region around CAL 83. The Hen emission peaks at 
the position of CAL 83, but it has a distinctly asymmetrical dis- 
tribution around the central star. We did not detect any He i emis- 
sion. 

We estimated an average value for the electron density n e in 
our FoV of -10 cm -3 which is consistent with the value found 
by RRM95. 

In addition to the flux maps, we also show flux ratio maps 
which characterise zones with different ionisation degrees in our 
FoV. The flux ratio values of four interesting positions in our 
FoV were used to place CAL 83 in the diagnostic diagrams of 
Veille u~& Osterbrock| ( 1987 ) which is used to distinguish be- 
tween AGNe and H n region-like objects (i.e. star forming galax- 
ies). This has not been done before for a SXS, as CAL 83 is the 
only know SXS surrounded by an ionisation nebula. We show 
that CAL 83 has characteristics of both AGNe and H n region- 
like objects and does not seem to be distinguishable from ULXs. 
We have also placed CAL 83, both observed and modelled, on a 
plot of [O m] , 15007/Hff vs He n ,14686/Hfi, developed especially 
for SXSs by |Rappaport et al.| {1994a] ). CAL 83, observed and 
modelled, differ slightly in the Heii/14686/H/3 ratio and more 
in the [O in] A5001/Hfi ratio, and is situated with the ULXs and 
starburst galaxies rather than with the Hn regions. 

Finally, we have presented a comparison between our obser- 
vations and model calculations for nebulae surrounding SXSs 
using the CLOUDY ionisation code. Improvement on the mod- 
els of Rappaport et al. (1994a) was achieved by utilising LMC 
abundances rather than Galactic. We found that none of the mod- 
els presented matches our observations completely. Keeping in 
mind that the observed nebula is not even close to being homoge- 
neous, as is assumed in the modelling, one could argue that the 
models are perhaps too simplistic. The modelling inconsisten- 
cies do not, however, affect the observationally oriented results 
presented in this paper. Perhaps the discrepancies are telling us 
something interesting and potentially important about the phys- 
ical processes that we do not yet understand. To get to the bot- 
tom of this, new observations covering the entire CAL 83 nebula 
should be compared to more detailed models which perhaps take 
into account a certain degree of inhomogeneity. 

Acknowledgements. The authors wish to thank the VIMOS support staff at 
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X 
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Model 5 


Model 6 


Observed 
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HI 
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3.10 
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HI 
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0.47 


0.46 


0.46 
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0.03 
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0.04 
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0.05 
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0.96 


6731 


[S II] 
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1.21 


1.74 


2.67 


0.48 


0.19 


0,65 



Notes. The Hj3 -normalised dereddened intensities of the emission lines taken from Table 3 compared to the cloudy model intensities. Wavelengths 
(A) are given in A and upper limits are indicated with a <. 
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